Microcanonical Foundation for Systems 
with Power-Law Distributions 
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Starting from microcanonical basis with the principle of equal a priori probability, it 
is found that, besides ordinary Boltzmann-Gibbs theory with the exponential 
distribution, a theory describing systems with power-law distributions can also be 
derived. 
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Systems exhibiting (quasi) power-law behavior in their probability distributions are 
quite ubiquitous in nature [1,2]. It is a challenging problem to understand the properties 
of such systems based on principles of statistical mechanics, because these systems can 
hardly be described by ordinary Boltzmann-Gibbs canonical ensemble theory, whose 
distribution is of the exponential form. In this Letter, we show that it is actually possible 
to derive non-Boltzmann-Gibbs theory from the microcanonical basis with the principle 
of equal a priori probability. The resulting distribution is found to exhibit the required 
power law. 

We first recall the Gibbs theorem [3,4], which states that a small subsystem of a 
microcanonical ensemble with large degrees of freedom is uniquely characterized by 
the standard exponential distribution. Historically this theorem has repeatedly been 
proved in various ways, such as the method of counting, the method of steepest 
descents, and the application of the central limit theorem. Starting from the 
microcanonical basis with principle of equal a priori probability, one obtains the 
exponential distribution for the canonical ensemble. Therefore, if the Gibbs theorem 
were universal, then any of the equilibrium theories other than Boltzmann-Gibbs theory 
could not exist and consequently power-law distributions would be excluded, as long as 
microcanonical ensemble theory is the basis. Here, we show that a route to canonical 
ensemble theory for systems with power-law distributions is actually allowed within the 
microcanonical ensemble theory. 

To exhibit this route, we begin with the standard discussion [5] about the Gibbs 
theorem by considering a classical system s and take its N replicas s v s 2 ,L , s N . The 
collection S = { s a} a=1 2 L n 1S re f erre d to as a supersystem. Let A a be a physical 
quantity associated with the system s a . It is a statistical random variable and its value is 
denoted by a (m a ) , where m a labels the allowed configurations of s a . The quantity of 
physical interest is the average of {.A a } _ l2h N over me supersystem: (1/N)^ a _j A a . 

According to microcanonical ensemble theory, the probabilities of finding S in the 



configurations in which the values of the average quantity lies around a certain value a , 
i.e., 
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are all equal. £ is assumed to be 



o(n- 1 - 8 ).. 
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where 8 > 0. The equiprobability P [m v m 2 , L , m N ) associated with this condition is 
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where 9(x) stands for the Heaviside unit-step function. To shift from microcanonical 
ensemble theory to canonical ensemble theory, we fix the objective system and 
eliminate the other ones. The probability of finding the objective system, say s v in the 
configuration m l = m is given by 



p(m)= P {m,m 2 ,L ,m N ), 



(4) 



m 2 , L ,m N 



which characterizes the canonical ensemble. 

To prove the Gibbs theorem, the following integral representation of the step 
function is employed: 



P + <j}X 

9(x)= f , 
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where (3 is an arbitrary positive constant. Then, in the large- N limit, one applies the 
method of steepest descents to evaluate the integration over (p. Contextually, it is clear 
that the exponential distribution in Boltzmann-Gibbs canonical ensemble theory has its 
origin in this integral representation of 0(x) using the exponential function. 

To examine the possibility of obtaining the power-law-type distribution, we consider 
the "^-exponential function", which is given by 



In the limit q — > 1, e q {x) converges to the ordinary exponential function. An important 

point is that even if the exponential function in the integrand in eq. (5) is replaced by 
the ^-exponential function, still the equality 



e q (x) = [l + (l-q) X ] 



(6) 



where q is a real parameter satisfying the condition 



q>\. 



(7) 
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as long as (5 is taken to satisfy 



l-{q-l)p XrmK >0. 



(9) 



The simplest way to ascertain this fact may be to use the following integral 
representation: 



(<px) = — — ^Adtt qX exp{-[l + (l-#)0x]*}, 



(10) 



where Hs) is the gamma function. The condition in eq. (9) will be discussed later in 
the context of the steepest descent approximation. 

Equation (8) may be understood from the point of view that the step function is of 
discrete topology and therefore can remain invariant under continuous deformation of 
the exponential function in the integrand in eq. (5). In this respect, it is clear that the use 
of the ^-exponential function is nothing but one particular choice of deformation. 

Using eq. (8), we have, in the leading order in N , 
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To obtain this equation, first we factorize the ^-exponential function as e q (^)[£ - M]) ~ 
e q {(j)£) e q {-(j)M) with M = (1/ N)^ a(m a )-d and then further factorize e q {-(j>M) 
into the product over the supersystem. Since | M \ < £ ~ o{n~ 1 ~ s } due to eqs. (1) and (2), 

this manipulation will be justified in the subsequent discussion of the method of 

steepest descents in the large- ./V limit. 

Noting that 0{£ - 1 M |) = 6(£ — M) — 9(—£ — M) and changing the integration 

variable as <fi — > Nip, we obtain 
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where sinn^x) = \e q (x) - jl and 
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Here, we examine the condition in eq. (9). In the present context, it is written as 



l-(<7-l)j3*Ar|±£-M| max >0. (14) 

The rectangular distribution function we are considering has a very narrow support. In 
fact, | ± £ - M | max is a quantity of o{N~ l ~ s } with 8 > . Therefore, in the large- ./V limit, 

ft* can be an arbitrary positive constant. 

Thus, working out to the leading order in N, we can express the probability as 
follows: 



P q ( m ) = p ( m > m 2> L > L m N ) 
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where W is the number of possible configurations satisfying eq. (1) and 

Z>) = L^(-0[a(m)-a]). (16) 

m 

Clearly, 

P 'T sinh (A^e) r ~ , ,i 
W= j # ^-^expiVInZ 9 (0)], (17) 



since p g (m) is normalized. Now, using the real part f3* of (j), the steepest descent 
condition reads 



— f = 0, (18) 
dp 



which leads to 



Pq{m) = iff) e ^~^ a{m) ~ ^' (19) 
a = Y, P q (m)a(m), (20) 



simultaneously, where P (m) is given by 



These results follow from eq. (6) and the relation de q {x)j dx = [e g (x)Y . It is evident 
that in the limit q — > 1 + all the discussions become reduced to the ordinary ones in 
Boltzmann-Gibbs theory with the familiar canonical distribtuion of the exponential 
form. On the other hand, the distribution function in eq. (19) is seen to asymptotically 
exhibit the power-law behavior 



[a{m)\ 



as desired. 



In the field of thermodynamics of chaotic systems [6], P q {m) in eq. (21) is referred to 

as the escort distribution, which is also a probability distribution associated with the 
original distribution p q (m). The steepest descent condition yields the fact that the 
arithmetic mean of 2 L N coincides with the generalized expectation value with 

respect to the escort distribution as in eq. (20). 

In conclusion, we have shown that not only ordinary Boltzmann-Gibbs canonical 
ensemble theory but also a theory for systems with power-law distributions can be 
obtained from microcanonical basis with the principle of equal a priori probability. It is 
worth pointing out that the structure in eqs. (19)-(21) is the characteristic of 
nonextensive statistical mechanics [7]. It is known that generalized canonical ensemble 
theory derived from maximum entropy principle based on the Tsallis entropy [8] with 
eq. (20) as the constraint gives rise to the distribution in eq. (20). However, it is 
essential to recall that here we made no initial assumptions on the definition of 
statistical expectation value and the form of the entropy. 
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